Degenerate Landau-Zener model: Exact analytical solution 
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The exact analytical solution of the degenerate Landau-Zener model, wherein two bands of de- 
generate energies cross in time, is presented. The solution is derived by using the Morris-Shore 
transformation, which reduces the fully coupled system to a set of independent nondegenerate two- 
state systems and a set of decoupled states. Due to the divergence of the phase of the off-diagonal 
element of the propagator in the original Landau-Zener model, not all transition probabilites exist 
for infinite time duration. In general, apart from some special cases, only the transition probabilities 
between states within the same degenerate set exist, but not between states of different sets. An 
illustration is presented for the transition between the magnetic sublevels of two atomic levels with 
total angular momenta J = 2 and 1. 

PACS numbers: 32.80.Bx, 33.80.Be, 03.65.Ge, 34.70.-|-e 



I. INTRODUCTION 

The Landau-Zener (LZ) model jl| is the most popu- 
lar tool for estimating the transition probability between 
two states whose energies cross in time, a situation which 
can occur in virtually every area of quantum physics. The 
LZ Hamiltonian involves the simplest nontrivial time de- 
pendence: a constant interaction and linearly changing 
energies. Nevertheless, due to some mathematical sub- 
tleties, when applied to real physical systems with more 
sophisticated time dependences the LZ model often pro- 
vides more accurate results than anticipated. 

The popularity of the LZ model, largely ensuing from 
the extreme simplicity of the transition probability, has 
stimulated numerous extensions to multiple levels. There 
are two main types of generalizations: single-crossing 
bow-tie models and multiple-crossings grid models. 

In the bow-tie models all energies cross at the same 
instant of time. Analytic bow-tie solutions have been 
found for three [31 and N states 0, i, i]. Examples of 
such systems occur, for instance, in a rf-pulse controlled 
Bose-Einstein condensate output coupler [1, 0|, and in 
the coupling pattern of Rydberg sublevels in a magnetic 
field 4i]. An extension, where one of the levels is split 
into two parallel levels, has been solved by Demkov and 
Ostrovsky [s*]. 

In the multiple-crossings models, a set of Na paral- 
lel equidistant linear energies cross another set of N\j 
such energies, thus forming a grid of crossings (Demkov- 
Ostrovsky model) [1 [IS [13 . For iV^ = 1 (or Na = 1) the 
Dcmkov-Ostrovsky model reduces to the earlier Demkov- 
Osherov model [ij, [3] ■ The special case when Na and 
Nb are infinite (so that the grid of crossings is fully peri- 
odic) has also been solved [l4!| . In the most general case 
of an arbitrary linear Hamiltonian, H(i) = D + Ci, where 
C is diagonal, the general solution has not been derived 
yet, but exact results for some surviva,l probabilities have 
been conjectured [5] and derived [isl. Ha. [TtI [Tsj . 

A variety of physical systems provide examples of mul- 
tiple level crossings. Among them we mention ladder 



climbing of atomic and molecular states by chirped laser 
pulses harpoon model for reactive scattering [2l|, 

and optical shielding in cold atomic collisions [2^ . 

A general feature of all soluble multilevel crossing mod- 
els is that the transition probabilities Pm^n between 
states 1pm and ipn are given by very simple expressions, 
as in the original LZ model, although the derivations are 
not trivial. In the grid models, in particular, the exact 
probabilities Pm^n have the same form (products of LZ 
probabilities for transition or no-transition applied at the 
relevant crossings) as what would be obtained by naive 
multiplication of LZ probabilities while moving across the 
grid of crossings from ipm to 'i/'n, without accounting for 
phases and interferences. For instance, the counterintu- 
itive transitions, for which the level crossings appear in a 
"wrong" order in time, are forbidden at infinite times. It 
has been shown, though, that the probability for coun- 
terintuitive transitions is nonzero for finite interaction 
duration [2^ or for piecewise-linear sloped potential [13] ■ 

An interesting feature of the existing multistate LZ so- 
lutions is that the respective derivations (usually using 
Laplace transforms and contour integration) all fail in 
the limit of degenerate levels and the assumption of non- 
degeneracy is essential. Effects of level degeneracies in 
the Demkov-Osherov model have been studied by reduc- 
ing the multistate dynamics to that of a single nonde- 

Pierate two-state system and several decoupled states 
, d^. Effects of quasi-degeneracies have been de- 
scribed by treating a nondegenerate system with small 
energy gaps as a perturbed degenerate system Q . 

In this paper, we derive the exact analytical solution 
for two crossing degenerate levels a and 6, of arbitrary 
degeneracies Na and Nb, which we shall refer to as the 
degenerate LZ model. Our model can therefore be consid- 
ered as an extension of the standard nondegenerate two- 
state LZ model to two degenerate levels. It also general- 
izes the solutions by Yurovsky and Ben-Reuven [25j and 
by Kyoseva and Vitanov [2^, which assume one degen- 
erate and one nondegenerate level. This model can also 
be viewed as the unsolved limiting case of the Demkov- 
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Ostrovsky model [1, [13, [lH for vanishing level spacing. 
Finally, this model represents the unsolved limiting case 
of the bow-tie models [1, [3, S 11] when all energies, that 
cross at the same time, coalesce into only two different 
slopes. 

Our method of solution is drastically different, and 
much simpler than those used in the nondegenerate mul- 
tistate LZ models. We make use of the powerful Morris- 
Shore (MS) transformation, which reduces the dynamics 
of two sets of degenerate states into that of a collection 
of independent nondegenerate two-state systems and de- 
coupled (dark) states. Each of the independent two-state 
systems represents a standard, nondegenerate LZ prob- 
lem, whereas the decoupled states do not evolve. Hence 
the solution of the degenerate LZ problem is equivalent 
to a collection of two-state LZ solutions. However, the 
situation is not so trivial because the different LZ so- 
lutions interfere and produce interesting features in the 
probabilities in the original basis. In particular, it turns 
out that not all transition probabilities are defined, as far 
as an infinite interaction duration is concerned. 

Among the numerous possible physical realizations of 
the degenerate LZ model, we point out the degenerate 
two-level system formed between two atomic levels of an- 
gular momenta J a and Jf, — J a or ± 1, driven by lin- 
early chirped laser fields of arbitrary polarizations. In 
the absence of magnetic field such a system represents 
exactly a degenerate LZ model. 

This paper is organized as follows. We define the 
problem in Sec. |TT] and the propagator is derived in 
Sec. mil in the general case. A special example for a 
Ja = 2 ^ Jb = 1 transition is considered in Sec. IIVI The 
conclusions are summarized in Sec. |Vl 
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FIG. 1; (Color online) The Morris- Shore transformation; a 
multistate system consisting of two coupled sets of degenerate 
levels is decomposed into a set of independent nondegenerate 
two-state systems and a set of decoupled states. 



We adopt a state ordering wherein the Na sublevels of 
the a level are placed first, followed by the N\j sublevels of 
the b level. In the rotating-wave approximation (RWA) 
the Schrodinger equation of the system reads [22] 



II. DEFINITION OF THE DEGENERATE 
LANDAU-ZENER MODEL 



ih-c{t)^n{t)C{t), 



(2) 



We consider a quantum system with Na degenerate 
states {|i/'m)}m'Li ™ the (lower) a set and TVf, states 

{\4'Na+n)}'^'Li i'^ the (upper) h set, as displayed in Fig. 
[T] (top). Without loss of generality we assume that 
Na ^ Nb. Each of the a states |'0rn) is coupled to each 
of the b states \^n) by a constant couphng Omn, and all 
couplings can be different. The a states are not coupled 
to each other directly, neither are the b states. All fields 
are off resonance by the same detuning A(i), which is 
assumed to be linear in time, with a rate C {chirp in 
coherent atomic excitation [23, HI]), 



^mn = const, 

A(i) = Ct, 



(la) 
(lb) 



For iVfc = 1 the present model reduces to the iV-pod 
model solved by Kyoseva and Vitanov [2^ , and for Na ~ 
A^b = 1 to the nondegenerate original LZ model. The 
present model is therefore a generalization of the A^-pod 
model to a degenerate b level. 



where the elements of the [Na + A"b)-dimensional vector 
C{t) are the probability amplitudes of the states. The 
adopted state ordering allows us to write the RWA Hamil- 
tonian as a block matrix, 



H(t) = 



V 

Vt D(t) 



(3) 



Here is the A^a-dimensional square zero matrix, in 
which the zero off-diagonal elements indicate the absence 
of couplings between the a states, while the zero diago- 
nal elements show that the a states have the same energy, 
which is taken as the zero of the energy scale. The ma- 
trix D(i) is an A'b-dimensional square diagonal matrix, 
with A(t) on the diagonal, D(t) = /^.{t)!^^- The absence 
of off-diagonal elements in D reflects again the absence 
of couplings between the b states, while the diagonal el- 
ements A stand for the common energy of all b states. 

In Eq. ([3]) V is an [Na x A'ft)-dimensional interaction 
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matrix with constant elements, 



flu fli2 



n 



^[\n,),\n2),...,\^Nj], 



(4) 



where |ri„) (n = 1, 2, . . . , Nh) are A^a-dimensional vectors 
comprising the interactions of the nth state of the b set 
with aU states of the a set, 



^2n 



n 



(n= 1,2, 



(5) 



III. EXACT ANALYTIC SOLUTION OF THE 
DEGENERATE LANDAU-ZENER MODEL 

A. Morris-Shore transformation 



decomposition of H into a set of independent two-state 
systems requires that, after removing the null rows, V re- 
duces to a iVfj-dimcnsional diagonal matrix; let us denote 
its diagonal elements by A„ {n — 1,2, ... , Nb). 

It follows from Eq. ^ that 



VV^ = AVV^At, 
V^V = BVtVBt. 



(10a) 
(10b) 



Hence A and B arc defined by the condition that they 
diagonalize VV^ and V^V, respectively. It is important 
to note that the square matrices VV^ and V^V have 
different dimensions, Na and Ni,, respectively. Because 
all elements of V are constant, A and B are also constant. 
It is straightforward to show that the Ni, eigenvalues of 
V^V are all non-negative; according to Eqs. ([9]) and 
(fTO| they are (n = 1, 2, . . . , Nf,). The matrix VV^ 
has the same eigenvalues and additional Nd = Na — Ni, 
zero eigenvalues. 

The MS Hamiltonian ([8| has the explicit form 



We shall solve the degenerate LZ problem by using 
the Morris-Shore (MS) transformation [29]. Morris and 
Shore have shown that any degenerate two-level system, 
in which all couplings share the same time dependence 
(constant in our case) and the same detuning (linear 
here), can be reduced with a constant unitary transfor- 
mation S to an equivalent system comprising only inde- 
pendent two-state systems and uncoupled (dark) states, 
as shown in Fig. [1] This transformation reads 



[A)=J2^*M, (6) 



where the tildas denote the MS basis hereafter. The con- 
stant transformation matrix S can be represented in the 
block-matrix form 



A O 
O B 



(7) 



where A is a unitary A^a-dimensional square matrix and 
B is a unitary iVb-dimensional square matrix, AA^ = 
At A = Ijv^ and BB^ = B^B = Ia,,. The constant 
matrices A and B mix only sublevels of a given level: A 
mixes the a sublevels and B mixes the b sublevels. The 
transformed MS Hamiltonian has the form 



H(t) = SH(i)St = 



O V 

Vt D(i) 



where 



V = AVBt. 



(8) 



(9) 



The Na X Nb matrix V has Nd = Na - Nb null rows 
{Na ^ Nb), which correspond to decoupled states. The 



H 













Ai 
A2 









Ai 
A, 




A 
A 















(11) 



The structure of H shows that in the MS basis the dy- 
namics is decomposed into sets of Nd decoupled single 
states, and Nb independent two-state systems ^ 
iV'n) {n = 1,2, . . . , Nb), each composed of an a state 
and a b state \tpn)i ^^'^ driven by the Hamiltonians, 



Unit) 



A„ 
A„ A(t) 



(n= 1,2, 



.Nh 



(12) 



These two-state Hamiltonians have the same detuning 
A{t) but different couplings A„. Each of the new a states 
iV'^i) is the eigenstate of VV^ corresponding to the eigen- 
value A^, whereas each of the new b states \ip'^) is the 
eigenstate of V^V, corresponding to the same eigenvalue 
A^. The square root of this common eigenvalue, A„, rep- 
resents the couphng between \tp!^) and \tpt.)- The Nd 
zero eigenvalues of VV^ correspond to decoupled (dark) 
states in the a set (since we assume that Na^ Nb, dark 
states, if any, are in the a set). The dark states are de- 
coupled from the dynamical evolution because they are 
driven by one-dimensional null Hamiltonians. 
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B. Solution to the degenerate LZ problem 

1. The MS transformation 

The MS decomposition allows us to reduce the degen- 
erate two-level LZ problem to a set of nondegenerate two- 
state LZ problems, wherein the detuning is unchanged 
and given by Eq. (fTb|) while the couplings A„, defined as 
the square roots of the eigenvalues of V^V, are combina- 
tions of the initial couplings between the a and h states. 

From the vector form (|4]) of V we obtain 



n=l 



vtv = 



{^2\^N,) 



(13a) 



(13b) 



Note that V^V is the Gram matrix for the set of vectors 
{\^n)}n=i- Thus if all these vectors are linearly indepen- 
dent then det V^V ^ and all eigenvalues of V^V are 
nonzero (30l |: however, this assumption is unnecessary. 

We assume that we can find the eigenvalues A^^ (n = 
1,2,..., Nb) of the matrices (|13ap and (jlSbp . and the 
corresponding orthonormalized eigenvectors: the Ni, cou- 
pled eigenstates |V'JJ) of VV''' and IV'Ji) of V^V, and the 
Nd decoupled eigenstates IV^^) of VV^ We use these 
eigenstates to construct the transformation matrices as 



B 



(14) 



Then according to the general theory the transformed 
interaction matrix ([9|) in the MS basis takes the form 
(fTTjl , where the positions of the Nd zero eigenvalues and 
the Nh eigenvalues A„ are determined by the ordering of 
the eigenstates in the transformation matrices (HD). 



2. The MS propagators 

Because the dark states are decoupled and have zero 
energies, their propagator is the unit matrix Ijv^. 

The propagator for each of the two-state MS Hamilto- 
nians is the LZ propagator for the respective cou- 
pling Xn, 



U„ = 



1 



6 = Ait)dt = -{ri- r, 



(15a) 
(15b) 



The Cayley-Klein parameters are [3l| 



AK^(r/e-'"/'^)A.^-i(r.e3"/4) 



-AK^(r/e3^^/4)A.^-i(r.e-*'^/4) 



(16a) 



t/4 



-fAK^(r/e3"/*)A.^(T.e-"/4) 



(16b) 



/3 = 



where k = f^o/vC, r — WC, and D^{z) is the parabolic- 
cylinder function. Ti = tiVC and Tf — tf^/C are the 
scaled initial and final times, respectively. In the original 
LZ model, — > —oo and Tf —f oo, and the Cayley-Klein 
parameters read [3l| 



a„ 

I3„ 



o-2TrA„ 



with 



:LZ 



c 

^2 



-A„ In ( 



iLZ 



- + argr(l 



(17a) 
(17b) 



(18a) 

(18b) 
(18c) 



Hence the phase diverges, which is a result of the 
unphysical assumption of an infinitely long interaction 
duration. This divergence is unimportant in the original 
LZ model because the transition probability. 



P. 



1 



-27rA„ 



(19) 



is well defined. Hence the final populations are well de- 
fined if the system starts in one of the two states, which 
is usually the case. However, when the system starts in a 
superposition of states, this divergence does not allow to 
calculate the populations, even in the original LZ model. 
We shall show below that in the degenerate LZ model 
this divergence does not allow for definite values of some 
populations even when the system starts in a single state. 

There are two divergent terms in the phase (|18b[) : 
polynomial and logarithmic, with different origins and 
different implications. The term \{Tf + Tj) is unimpor- 
tant in the present context because it derives from the 
chosen Schrodinger representation ([3]); in the interaction 
representation (when the detunings turn into phase fac- 
tors of the couplings) it disappears. Moreover, this term 
is the same for all /3„ and factors out of the probabilities 
(see below). The term iA„ ln(T,j^T|), however, depends 
on /3„; it arises from the nonvanishing coupling and the 
rather slow divergence of the detuning. These logarith- 
mic terms cannot be factored out, unless the MS cou- 
plings A„ coincide or vanish by accident, and appear in 
some transition probabilities, as we shall see below. 
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The propagator in the original basis 



By taking into account the LZ propagators ([15)) for the 
Ni, two-state MS systems, the ordering of the states, and 
the MS Hamiltonian pT|) . the fuU propagator in the MS 
basis can be written as 



U = 
















ai 


• 


• 


-PI 


• • 








a2 ■ 


• 





-(3*2 •• 








• 


• OiNt 





• • 






• 


■ 


a* 


• • 








(32 ■ 


• 





a*2 ■■ 








• 







• • 





By using the completeness relation 



(20) 
(21) 



it is straightforward to show that the propagator in the 
original basis U = StUS reads 



U 



i+E^:iK-i)IC)(CI 



(22) 

Note that the propagator does not depend on the decou- 
pled states {k = 1, 2, ... , Nd), which are excluded 
by using Eq. ([2T|) . This has to be expected because, 
owing to their degeneracy, the choice of the decoupled 
states is not unique: any superposition of them is also 
a zero-eigenvalue eigenstate of VV^. Because the dy- 
namics in the original basis must not depend on such 
arbitrariness, the propagator U must not depend on the 
decoupled states at all. 



4. Transition probabilities 

If the system starts in an arbitrary state of the a 
set then Eq. ([22]) gives for the matrix elements Ufi = 

(tpflUl^pi) the expressions 



E„=i ("n - 1) a/nflm (/ ^ a Set), 

X/n = l (3nbfnO-in if e b set). 



(23) 

where akn and b^n denote the components of the MS 
states \tp'^) and respectively, 

IV'n) = [ain,a2n,---,ajv„n]'^, (24a) 

\i^n) = [bln,b2n,..-,bN,nf . (24b) 

If the initial state belongs to the b set, we have 

TT _ / - /3>/«&m (/easet), , 

E:U<bf^bL ifebset). ^''^ 



In both cases, the transition probability from state 
to state IV'/) is 



P^^f - \U 



(26) 



Equations ([23)) - (l26|) reveal several important features 
of the degenerate LZ model. 

(i) The transition probability Pi^f is always well de- 
fined if the initial and final states belong to the same set 
of states (a or b) because then Pi^f involves only the 
Cayley- Klein parameters a„ , which are real and positive 
and hence do not have divergent phases, see Eq. ()17a|) . 

(ii) When the initial and final states belong to diff'erent 
sets, \ipi) to the a set and j?/"/) to the b set, or vice versa, 
the transition probability P-i^f is well defined only if the 
corresponding sums in Eq. (|23|) or ([25]) reduce to one 
term (because some of the a and b coefficients may vanish 
accidentally) or if the phases of all participating /3„'s are 
the same. The latter may only happen accidentally if all 
MS couplings A„ are equal: then the phases factor out 
and cancel in the transition probability. 

(iii) Baring accidental cases discussed in the previous 
point, the transition probabilities between states from 
different sets are not defined due to the divergence of the 
phases of the Cay ley-Klein LZ parameters 



5. Summary 



In summary, Eq. (1221) gives the propagator for the de- 
generate LZ model. The transition probabilities can be 
calculated from Eqs. ([23l) -([26 |) . which require the knowl- 
edge of the coupled MS states \'tp'^) of the a set and \^p^) of 
the b set. The former are the eigenstates of VV^ and the 
latter are the eigenstates of V^V. The knowledge of the 
decoupled zero-eigenvalue states of the a set is not 
necessary for the calculation of the propagator. Not all 
transition probabilities are defined for infinite time dura- 
tion because of the divergent phases of the Cayley-Klein 
parameters /?„. For any finite initial and final times, 
though, all transition probabilities are well defined. 



6. An alternative: the Allen-Eberly-Hioe model 

In a real physical situation with degenerate levels, a 
more realistic alternative to the LZ model is the lesser 
known AUen-Ebcrly-Hioc model [H, [s^l 



VLq sech(t/r), 
Btanh(t/r). 



(27a) 
(27b) 



Here the coupling VL{t) is a bell-shaped pulse, with a char- 
acteristic width T . The detuning crosses resonance at 
time t — Q and does not diverge at infinity but tends to 
the finite values ±5. The Cayley-Klein parameters for 
this model, including their phases, are well defined. 
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7 = 2 



I 



Morris-Shore 
transformation 



\W 



FIG. 2: (Color online) The Ja = 2 J;, = 1 transition in the 
original basis (top) and in the MS basis (bottom). With only 
circularly polarized fields the full eight-state system decouples 
into a five-state M system and a three-state A system. A 
linearly polarized field would couple the M and A systems. 
The Morris-Shore transformation turns the M system into a 
pair of two independent nondegenerate two-state systems and 
a decoupled state (bottom left), and the A system into a two- 
state system and a decoupled state (bottom right). 



IV. EXAMPLES 

A. Ja ~ 2 ^ Ji, = 1 transition 

1. General case 

We illustrate the above results with a specific example: 
transition between two atomic levels with total angular 
momenta Ja — 2 and Jf, = 1 in the field of two circu- 
larly polarized (right and left a~) chirped-frequency 
laser fields with linear chirp and steady amplitudes. In 
the absence of magnetic fields, the 5 magnetic sublevels 
of the Ja level are degenerate and so are the 3 magnetic 
sublevels of the Jt level, as shown in Fig. [2l This system 
therefore represents a physical realization of the degen- 
erate LZ model with Na = 5 and Ni, = 3. When only a'^ 
and polarized fields are present the eight-state system 
decouples into a five-state M system, which is composed 
of the sublevels with Ma = —2,0,2 and Mb = —1,1, 
and a three-state A system comprising the sublevels with 
Ma = — 1, 1 and Mb = [1^. If there is also a linearly 
(tt) polarized field then the M and A systems couple and 
all eight states will be involved in the dynamics. 

The two and a~ polarized fields can be produced 
by a single elliptically polarized field; then the amplitude 
ratio and the relative phase of the and a~ fields can be 
controlled, respectively, by the ellipticity and the rotation 



TABLE I: Coefficients of the MS basis states (|3TJ for the 
Ja = 2 <-> Ji, = 1 transition. The values for arbitrary ellipti- 
cal polarization e are in the second column, and those for a 
linear polarization £ = in the third column. The relevant 
normalization coefficient i/ for elliptical polarization is listed 
after each group of coefficients. 



d'-2 
do 



a_ 



2,1 



0-0,1 

4,1 
2 



la 



a-2,2 
4,2 



6-1,1 
bi,i 

6-1,2 
6l,2 



arbitrary e 



/^d(l 



(1 + e) 
4^2-e2^ 



-\via{l+e) (l-6e- Vl + 24e2) 
viae^J&{l - e^) 



\via (1 - e) (1 4- 6e - VI + 24e2) 
VI + 24£2 [(1 + ^1 + 24£2 + (lle^ ^ 1)] 



-\v2a{l+e) (1 - 6£ -f- VI + 24£2) 



\v2a (1 - £) (1 + 6£ + Vl+24£2) 
VI + 24£2 [(1 + £2) VI -I- 24£2 - (ll£2 - 1)] 



//6\/VrT24£2 + 5£ 
-/^6\/Vl + 24£2 + 5£ 



2V1 + 24£2 



£ = 



1 

2 



7i 



angle of the field. Moreover the (7+ and a fields will 
have automatically the same detuning. 

We shall only consider the M-system, because the A- 
system contains a non-degenerate upper state and can be 
treated with a simpler formalism [20| . 

The interaction matrix for the M system, with the 
Clebsch-Gordan coefficients accounted for, reads 1331 



V 



1 



y^n+e''>+ 
V6^-e" 



(28) 



and hence the matrices VV^ and V^V are 

vvt 



1 

10 



(29a) 



1 

10 



6nl + n^_ o+o_e'^ 



(29b) 



with 6 = 6^—9^ being the relative phase of the two fields. 
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The eigenvalues of VV^ are (n = 0, 1, 2), where 

Ao = 0, (30a) 



Al 2 = ^\ 



'7± Vl + 24£2 
20 '' 



(30b) 



with n = ^Jnl + fll and e = + ^l)- 

The eigenvalues of V^V are A^ and Aj. The eigenstates 
of VV^ are a decoupled state and two coupled states 
iV'i) and IV'2)' which are composed of a states, whereas 
the eigenstates of V^V are two new b states j33i] . 









m=-2,0,2 








m=-2,0,2 




m=-l,l 



ime/2 



(31a) 



The coefficients of these new MS basis states are given 
in Table U [33,] ; they are related to the coefficients in Eqs. 
(ESD-dSSl) as dm = <e-™^/2^ arar, = o!mn&-'""^ 1'^ , and 
&mri — &mnC~™^^^- By using thcsc coefficients and Eqs. 
([lil)-([TH]), (HSD-dlSl), and (1201), one can find the transition 
probability between any two states. 



2. Case of 



trigs 



We shall consider in some detail the special case $7+ = 
f2„; then e = and the coefficients in Table [J simplify. 
The MS couplings (|30b[) . the LZ factors, and the Cayley- 
Klein parameters reduce to 



Al 



ai 
a2 




(32a) 

(32b) 

(32c) 
(32d) 



where ^ = TrfJ^/lO/?. It is particularly significant that 
the coefficient ao,2 associated with state |'0o) vanishes 
accidentally, ao.2 = 0, see Table ID The implication is 
that the sums in Eqs. (|23p and (|25p . which involve 
coefficients, reduce to just single terms when state |^o) is 
involved. Consequently, all transition probabilities from 
and to state |'0o) are defined and the divergence of the 
phases 0„ does not show up here. 

The propagator in the original basis reads (for 9 — 0) 



U = 
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V 32 
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" (33) 

In the adiabatic limit ^ ^ 1 the matrix P = {-P/i}i,/=-2.o,2.-i,i with the transition probabilities Pi~,f — Pfi reads 
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The /3's which are left over in the propagator ((33|) have 
divergent phases, see Eqs. (|17bp and (Il8bj) . Because the 
respective couplings Ai and A2 [Eq. ([5^ ] are different, 
the logarithmic components in the phases of the /3's are 



different and therefore give rise to an interference term in 
the transition probability, which oscillates in time with 
a logarithmically increasing frequency. Hence the tran- 
sition probabilities with sums over different /3's do not 
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FIG. 3: (Color online) Time evolution of the populations in 
a five-state M system formed of the magnetic sublevels M — 
-2, 0, 2 of the J = 2 level and M = -1, 1 of the J = 1 level, 
for linear polarization (e = 0), in the cases when the system 
starts in (i) top: state | — 2) of the J = 2 level; (ii) middle: 
state |0) of the J = 2 level; (iii) bottom: state | — 1) of the 
J = 1 level. The arrows on the right indicate the asymptotic 
values at t — + oo, wherever applicable. The chirp rate C is 
used to define the time and frequency scales. The coupling is 
Q = 5C^''^, which implies that the adiabatic condition = 
2.5tv ^ 1) is fulfilled and the adiabatic solution (|34p applies. 
The initial time of the integration is ti = — 400C~^''^. 



have a limit at infinity. At any finite times, however, 
these probabihties are well defined. 

Figure [3] displays the time evolution of the populations 
of the five states in the ncar-adiabatic limit for linear 
polarization (e = 0) and for three different initial condi- 
tions. In the top frame the system starts in the J = 2 
state |'0-2)- As predicted by Eq. ([34|) the populations of 
the J = 2 states acquire definite values as t ^ oo, while 
the populations of the J = 1 states oscillate: the log- 
arithmic scale demonstrates that indeed, the oscillation 
phase diverges logarithmically. 

Figure [3] (middle frame) displays the time evolution 
of the populations when the system starts in the J — 2 
state IV'o)- As predicted by Eq. (p4)) the populations 
of all five states acquire definite values at infinity, that is 
all transition probabilities exist, because of the accidental 
vanishing of the coefficient 09,2, as discussed above. 

Figure [3] (bottom frame) displays the time evolution 



TABLE II: Morris-Shore couplings for transitions with Ja = J 
and Jj, = J — 1 or J for polarization e = 0. 



Ja = J and Jb = J — 1 



integer J 



2n(2J-n) 



J{2J + 1) 
(n = 0,l,..., J) 



half-integer J 



2n(2J-n) 



J(2J + 1) 
(n = 0,1,..., J -1/2) 



Ja — Jb 



J 



integer J 
-o 2n 

^27(7+1) 

(n = 0,l,..., J) 



(n: 



half-integer J 

-o + l 

\/'2J{J + 1) 
= 0,1,..., J -1/2) 



of the populations when the system starts in the J = 1 
state As predicted by Eq. the populations of 

the J = 1 states acquire definite values (zero) as t — > 00, 
while the populations of the J = 2 states oscillate, with 
a logarithmic divergence of the oscillation phase. The 
exception is the population of state \ipo), which exists 
because of the accidental vanishing of the coefficient ao,2- 



B. The case of arbitrary transition with Ja 

Jb — J — 1 or J 



J and 



For Ja = J and Jh = J — 1 with integer J, in the 
presence of right and left circularly polarized fields only, 
the full 4J-state system factorizes into two independent 
subsystems, like the M and A systems in Fig. [21 The 
larger, (2 J + l)-state system is formed of the magnetic 
sublevels Ma — —J,—J + 2,...,J of the Ja level and 
Mb = - J + 1, - J + 3, . . . , J - 1 of the Jb level. The 
smaller, (2 J — l)-state system is formed of the magnetic 
sublevels Ma = — J+1,—J + 3,...,J — 1 of the Ja level 
and Mb = - J + 2, - J + 4, . . . , J - 2 of the Jb level. For 
equally strong (t+ and a~ fields (e ~ 0) the MS couplings 
of the larger subsystem are given in Table [TTl The smaller 
subsystem has the same MS couplings, except for the 
largest one (with n = J). 

When J is half-integer the two independent subsystems 
are composed of similar sets of magnetic sublevels but 
with opposite signs of M . Because of this symmetry, the 
eigenvalues are exactly the same for both subsystems. 

For Ja = Jb = J the two subsystems are equivalent 
and they have the same eigenvalues, which are also listed 
in Table [TTl for both integer and half- integer J. 

The eigenstates (the MS states) are too cumbersome 
to be presented here, but they can easily be found for 
any particular J. 



V. CONCLUSIONS 

In this paper we have derived the solution of the 
time-dependent Schrodinger equation for the degenerate 
Landau-Zener model, which involves two crossing sets of 
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degenerate energies. The states in each set interact with 
the sublevels of the other set but there are no direct cou- 
phngs within the same set of states. A physical example 
is the transition between the magnetic sublevels of two 
levels with nonzero angular momenta induced by steady 
laser fields with linearly chirped frequencies. 

The solution uses the Morris-Shore transformation, 
which decomposes the original fully coupled system into 
a set of independent nondegenerate two-state LZ systems 
and a set of decoupled, dark states. Using the known two- 
state LZ propagators we use the inverse transformation 
to obtain the propagator in the original basis. 

Our results complement the Demkov-Ostrovsky model, 
which assumes two crossing bands of equidistant non- 
degenerate energies. Our results also complement the 
bow-tie models, which also exclude degeneracies. Our 
derivation is simpler than in these nondegenerate mod- 
els; however, the results are not so remarkably simple, 
as in these models, because of interferences between the 



different LZ propagators in the MS basis. More impor- 
tantly, we have found that not all transition probabilities 
exist for an infinite coupling duration, because this un- 
physical assumption gives rise to a divergent phase in the 
original nondegenerate LZ model. In the latter model the 
transition probability is not affected because this phase is 
cancelled. In the present degenerate LZ model, however, 
these divergent phases interfere and make some of the 
transition probabilities undefined in the limit of infinite 
times. As a rule, the transition probability between any 
two states within the same set always exists, but between 
two states from different sets can only exist by accident. 

Our results demonstrate that the LZ model should be 
used with care when multiple states are involved. In real 
physical situations the lesser known AUen-Eberly-Hioe 
model [28l.[33| can be a viable alternative, particularly in 
the presence of degeneracies, because it involves a pulse- 
shaped interaction, and hence no phase divergence. 
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